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The gauge-gravity duality can be used to relate connected multi-point graviton functions to con- 
nected multi-point correlation functions of the stress tensor of a strongly coupled fluid. Here, we 
show how to construct the connected graviton functions for a particular kinematic regime that 
is ideal for discriminating between different gravitational theories; in particular, between Einstein 
theory and its leading-order string theory correction. Our analysis begins with the one-particle 
irreducible graviton amplitudes in an anti-de Sitter black brane background. We show how these 
| can be used to calculate the connected graviton functions and demonstrate that the two types of 

■ amplitudes agree in some cases. It is then asserted on physical grounds that this agreement persists 

f"^ ■ in all cases for both Einstein gravity and its leading-order correction. This outcome implies that the 

£\J ' corresponding field-theory correlation functions can be read directly off from the bulk Lagrangian, 

. , just as can be done for the ratio of the shear viscosity to the entropy density. 

o. 

i. introduction 

t-h ; 

The gauge-gravity duality implies that a strongly coupled gauge field theory in four dimensions should have a dual 
description as a weakly coupled gravity theory in five-dimensional anti-de Sitter (5D AdS) space [H, 0. Because the 
gauge theory lives at the outer boundary of the AdS bulk spacetime, this is also a bulk-boundary or "holographic" 
O | correspondence [3J. 

String theory is a unitary and UV-complete theory. Consequently, any consistent effective description of the physics, 
on either of the side of the duality, must be unitary. This observation lead us to the following conclusion An 
effective theory describing small perturbations about the gravitational background solution cannot have more than 
t-H 1 two time derivatives in its linearized equation of motion. In other words, such an effective theory of gravity must 
\ either be Einstein's two-derivative theory, a member of the exclusive Lovelock class of higher-derivative theories Q 
or else be supplemented by boundary conditions that enforce this two-derivative limit. We will often use a term like 
"Lovelock" or "Gauss-Bonnet" (the four-derivative Lovelock theory) with the understanding that either of the latter 
two scenarios is realized. See Section 1.2 of [6j for further discussion. 

Lovelock gravity consists of an infinite series of terms that are arranged in order of increasing numbers of pairs of 
derivatives. But, given a sensible perturbative expansion, it is likely that the correction to the Einstein Lagrangian 
at leading order will be the most easily accessed by experiment (if at all). So that it makes sense to say that the 
gravitational dual is Einstein gravity plus a Gauss-Bonnet correction, or possibly just Einstein gravity. We will be 
proceeding with this as our premise. 

Our claims about the gravitational dual can be put to the test, as we now elaborate on. In a recent publication Q, 
we systematically studied gravitational perturbations about a black brane solution in 5D AdS space. We focused on 
a certain kinematic region of "high momentum" , which will be clarified early in the paper, but otherwise computed 
all of the physical, on-shell one-particle irreducible (1PI) graviton rt-point amplitudes for both Einstein gravity and 
its leading-order Gauss-Bonnet correction. Some of the more important findings are listed as follows: 

1) All 1PI n-point functions with n odd are trivially vanishing. 

2) All 1PI n-point functions involving the scalar graviton modes (or the sound channel |7|) are trivially vanishing. 

3) All 1PI n-point functions involving the vector graviton modes (or the shear channel [7|) are essentially featureless 
and theory independent. 

4) Except for the two-point function, any 1PI 2n-point function consisting strictly of tensor graviton modes depends 
on the scattering angles and in a theory-dependent way. 

5) For Einstein gravity, the (tensor) four- and higher-point functions all go as s, where s has its usual meaning as the 
Mandelstam center-of-mass variable. 

6) The contribution of the leading-order Gauss-Bonnet correction goes as s 2 for the four-point function and s(s + v) 
for higher-point functions, where v is an independent angle that can be interpreted as a "generalized Mandelstam" 
variable. 

In brief, what we have found is that both Einstein gravity and its leading-order Gauss-Bonnet correction exhibit 
a universal angular dependence in multi-point graviton functions (at least for a particular kinematic region). The 
two extra derivatives in Gauss-Bonnet gravity lead to a more complicated angular dependence, which would become 
increasingly more complex as the perturbative order of corrections is increased. But the added complexity only 
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becomes evident at the six-point (and higher) level and, even then, only for a certain class of gravitons. 

The black brane geometry provides a framework for probing the hydrodynamic properties of a strongly coupled fluid 
, such as the quark-gluon plasma . The gauge-gravity duality implies that there should be an analogous list 
of statements to the above list about the stress-tensor correlation functions of a strongly coupled gauge theory. This 
point was already mentioned in our previous work Q and has since been made concrete in a more recent article [l2j . 

An obstacle to directly applying the duality for this purpose is that the graviton amplitudes of @ are the 1PI 
amplitudes. What is rather needed to make contact with the gauge theory would be the connected graviton functions. 
The reason being that the 1PI functions on the field-theory side are not particularly useful because of complications 
that arise in separating these from the connected functions (see [HI, [l4| for a recent discussion). The best way to 
circumvent such issues is to have things already worked out in the bulk before applying the duality. 

The main objective of the current study is to fill in this gap on the bulk-to-boundary chain and make clear how 
the connected functions can be calculated from the 1PI ones. Working in the high-momentum kinematic region, we 
are able to show that the four-, six- and eight- point connected functions for Einstein gravity are the same as the 1PI 
ones. Similar agreement is established for the four- and six-point functions when the leading-order Gauss-Bonnet 
corrections are included. 

Following our presentation of explicit calculations, we will assert, using general covariance and additional general 
arguments, that this agreement between the connected and 1PI functions should persist for the higher-point amplitudes 
as well. This is an interesting outcome in its own right but is especially significant in the context of the gauge-gravity 
duality. 

What our findings make clear is that, at least for the high-momentum regime, the multi-point correlation functions 
of the gauge-theory stress tensor fall into universality classes [l2| in the same way that the ratio of the shear viscosity 
to the entropy density famously does And, just as one can read off this ratio of two-point correlators directly 

from the bulk Lagrangian [17Hl9j | , the same can be done for the higher-point correlation functions as well. Importantly, 
this provides a new means for probing the gravitational dual of a strongly coupled gauge theory. 

The rest of the paper proceeds as follows: The next section clarifies our choice of kinematic region and Section 3 
reviews some of the important results from our previous article @. In particular, we recall our formulations for the 
1PI graviton functions and discuss how these are related to correlation functions of the boundary theory. Section 4 
is where we begin to address the relation between the connected and 1PI amplitudes, starting with a preliminary 
discussion about the combinatorics of Feynman diagrams. In Section 5, we show by explicit calculation that the 
one-particle-reducible (1PR) contribution to the connected six-point function identically vanishes (the four-point case 
is trivial). This outcome is extended to all the reducible diagrams with a single internal line in Section 6. Then, 
in Section 7, we present our argument as to why these cancelations should universally persist. Section 8 considers 
disconnected graviton functions, which can also be useful for probing the gravitational dual. Section 9 ends the main 
text with a brief overview and Appendix A includes a calculation that is referred to in Section 7. 

Let us briefly mention coordinate conventions. We assume a 5D bulk spacetime and use x = (t, x, y, z) to denote 
the coordinates of the four-dimensional Poincare-invariant subspace of the brane. Without loss of generality, z is 
designated as the direction of graviton propagation along the brane (or along any other spacetime slice at constant 
radius). For indices, a, b, ■ ■ ■ — {t, x, y, z} . The radial coordinate is denoted by r and ranges from r = at the 
black brane horizon to r — > oo at the AdS boundary. 

II. THE HIGH-MOMENTUM REGIME 

Before continuing, we wish to clarify how the kinematic region of high momentum is defined. Here, the focus is on 
the bulk perspective; for the boundary point of view, see 

Essentially, the frequencies and momenta (w's and fc's) of the gravitons should be large enough so that these 
dominate over all radial derivatives, whether acting on the gravitons or the background. In practice, this amounts 
to requiring that any derivative acts on a graviton to produce an uj = —i\7 t = —idt or a k — — iV 2 = — id z . (39j 
This kinematic regime is most suitable for distinguishing between Einstein gravity and higher-derivative corrections, 
as it emphasizes the number of derivatives in an interaction vertex. 

The parameter e will serve as a dimensionless perturbative coefficient which measures the strength of the higher- 
derivative corrections to the Einstein Lagrangian, such that the Gauss-Bonnet correction is linear in e. [40( Then the 
number of derivatives is given by 2 1+ for a term of order e J . The restriction to high momenta is really the same as 
asking that the Mandelstam invariant s be larger than the AdS bulk curvature or s 3> 1/L 2 , where L is the AdS 
curvature length. 

As for an upper bound, the momenta should still be small enough to ensure that any higher-order Lovelock terms 
make parametrically small contributions to physical quantities. This is necessary for keeping the perturbative expan- 
sion under control and, more importantly, for consistency with treating string theory as an effective theory of gravity. 
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Now consider that each additional order in the Lovelock expansion brings with it another factor of e. For instance, 
a term with six derivatives can contribute one more factor of s but at the cost of one more factor of e. Hence, we 
require se <C 1/L 2 . And, because it follows on general grounds that e ~ I 2 ,/ Li 2 <C 1 (l p is the 5D Planck length and 
the right-hand side is a necessary requirement for the viability of the gauge-gravity duality) , the previous inequality 
means s <C 1/lp . That is, all momenta must be "sub-Planckian" . 

This last constraint is rather trivial, but there is still another one to consider. Because a black brane is dual to a 
fluid, we are working (at least implicitly) in the hydrodynamic regime, and so the bound s <C (jrT) 2 is also required. 
Here, T is the fluid temperature, which is equivalent to the Hawking temperature on the brane. We then end up with 
the high-momentum region being defined by the finite range 

1 < L^fs < irTL . (1) 

There is no danger of the right-hand inequality being in conflict with the bulk solution because the gauge-gravity 
duality can only be applied consistently in a black hole (or brane) background when the horizon is much bigger than 
the AdS scale L or, equivalently, much hotter than the surrounding AdS space [22|. Since ~ ttTL 2 , it follows 
that ttTL > 1 . 

An additional comment is in order: Let us recall the relation A = g YM N between gauge-theory parameters, where 
A is the 't Hooft coupling, gyu is the Yang-Mills coupling and N is the number of colors. The holographic dictionary 
applies when N — > oo and gyu ~^ , with A large but finite. Hence, we carry out our analysis at tree-level (loop 
diagrams are suppressed by inverse powers of N) and the leading-order modifications to Einstein gravity must then be 
at higher orders in A -1 . This is consistent with looking at higher-derivative extensions to the effective string-theory 
Lagrangian (23|. 



III. GRAVITON MULTI-POINT FUNCTIONS 

In this section, we recall some relevant parts of our earlier work [6| as well as set up the basic framework. 

The initial task is to calculate all the physical, 1PI on-shell graviton multi-point amplitudes for an (asymptotically) 
AdS theory of gravity. We assume that the background solution is a 5D black brane with a Poincare-invariant horizon 
because, as already mentioned, this setting is useful for learning about the fluid dynamics of the gauge-theory dual. 

As also mentioned, the effective theory describing small fluctuations about the background is constrained to have 
a two-derivative equation of motion. It is assumed that this effective theory limits to Einstein gravity in the IR and 
has a sensible perturbative expansion in terms of the number of pairs of derivatives. It follows that each term in the 
Lagrangian is suppressed by a factor of e relative to its predecessor with two derivatives fewer. Hence, to leading 
order, we need only consider the four-derivative correction to Einstein's theory; so that the theories of interest are 
Einstein gravity and Einstein gravity plus a Gauss-Bonnet extension. 

The premise is to expand the Einstein and Gauss-Bonnet Lagrangians in numbers of gravitons. A graviton h^ v 
represents a small perturbation of the metric from its background value, g^ v — > g^ + h^v . The n th -order term of the 
expansion can be identified with the graviton n-point function, which can be represented by a 1PI Feynman diagram. 
This becomes a Witten diagram Q in the limit that r — > oo . 

Let us introduce the following ansatz and labeling convention for the gravitons: hftfu = <f)(r) exp [iojjt — kjz] . We 
work in the radial gauge, for which h rr = h ar = for any choice of a. This choice leads to a decoupling of the 
gravitons into three distinct classes pj: tensors, vectors and scalars or (respectively) h% — {h xy } , hi — {h zx ,h tx } 
and ho = {htt, h zz , h zt , h x x + h y y } , with redundant polarizations omitted. 

As previously motivated, we work in the kinematic region of high momentum, meaning that only terms with the 
highest allowed power of uj and k are included. 

Given that the calculations are on-shell and in the region of high momentum, we are able to discard any n-point 
function involving scalar modes. This, in turn, means that graviton functions with odd values of n vanish, as general 
covariance dictates that vectors and tensors must come in pairs. 

The case against vectors modes is somewhat more subtle. A vector mode is analogous to a field-strength tensor [§[ , 
so it must be differentiated to be physical. Moreover, any non- vanishing n-point function that includes vector modes 
is limited to exactly two derivatives [6], both of which are necessarily constrained to act on the vectors. Because of 
this simplicity, this class of functions is of limited value to the current analysis and will subsequently be ignored. 

Thus, we are left with only the tensor modes. Later on, we will calculate connected functions for the tensors by 
convolving pairs of 1PI amplitudes, and so it is important to understand why vector and scalar modes can be ignored 
in the internal lines. Generally speaking, all types of modes can appear because, unlike those appearing on external 
lines, internal gravitons are not constrained to be on-shell. However, on-shell or off, the scalar modes are not physical 
(i.e., can be gauged away) in the absence of an external source. (See, e.g., |24j.) And, since the vector modes can be 
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identified with eletromagnetic gauge fields [8j , they must likewise be sourced to be physical. One might then ask how 
our results would differ when such sources are present, as could well be the case in a string-theory context. We will 
make an argument that, in the high-momentum regime and at least to order e, all internal lines - irrespective of the 
internal propagating modes - do not contribute. We will return to this matter at the end of Section 7. 

As explained in [6|, we have found the complete list of 1PI multi-point functions for the tensor modes at arbitrary 
r. But, for current purposes, it is most useful to focus on the limit of the AdS outer boundary; this being the location 
of the gauge-theory dual. In this case, 
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where E and GB respectively denote Einstein and Gauss-Bonnet corrected gravity, the 5D Newton's constant has 
been fixed according to 1Q-kG§ — 1 , and the kinematic variables s and v are made explicit below. We have also 
employed the boundary limit of the AdS brane metric, lim ds 2 = — j^dt 2 + ^dr 2 + \dx 2 + dy 2 + dz 2 ] . 



More generally, for 2n = 4, 6, 8 . 
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with the understanding that v = for 2n — 4 . 

The external gravitons are symmetrized in our expressions, with the center-of-mass variable s and the generalized 
Mandelstam variable v accounting for the symmetrization: 



^ Atl Alt 

2n(2n-l) J2J2 k i k ^ 

y 1 i=i 3=1 



(10) 



2n 2n 2n 



2n 

E 



2n(2n-l)(2n-2)(2n-3) ^ ^ ^ 

«2#il i3#il,2 24 7^*1,2,3 J#H,2,3,4 



2» 

2J ^ii ^> 



(11) 



The Gauss-Bonnet corrections depend only on the Riemann-tensor-squared term, as the other four-derivative terms 
(Ricci-tensor-squared and Ricci-scalar-squared) can be transformed away with a suitable choice of field redefinitions 
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[a, [2fl [26| ■ Also, to arrive at the above expressions, we employed the transverse/traceless gauge for the tensors (which 
is compatible with our previous choice of radial gauge). 

Let us briefly explain the above outcomes from a Feynman-diagram perspective: With knowledge that these are 
really the IPI graviton amplitudes, we are essentially computing a functional derivative of the form (see, e.g., [271 ]) 



</' 2 



H) 



2n+l 



f2n 



in 



J V[h] e 



i^gC[g,V.h]+iJh 



(12) 



where an h should be regarded as a tensor, J is an external source, there is an implied spacetime integral in the 
exponent and, since tensor modes come in pairs, the Lagrangian density in the exponent can be expanded out as 



h 2 h 4 h 6 

<J=g£\g,V,h] = C [g] + -^C 2 [g,V] + -^U[g,V] + -^C 6 [g,W] 

Now, taking the variations and then setting J = 0, we have 

JV[h] h ^V^c l9 ,vM 
[ ' JV[h] e V=9£b,v,h] • 

A Taylor expansion of the exponentials then leads to the outcome 

1 JV[h] h 2n (h 2n C 2n [g,V]) e* h * c *bW 
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(15) 



with the ellipsis indicating contributions that are not IPI. [4l| 

To proceed, one contracts together the pair of monomials (each of order 2n) within the upper integral of Eq. (|15l) . 
This can be done in (2n)! different ways, (42j | with this number canceling the graviton symmetrization factor in the 
denominator. After the standard process of amputating the external lines, the result of the contraction is simply 
£-2n{g, V). Equations (O-© are these quantities in momentum space. 

An important feature of these multi-point functions is how they depend on the scattering angles of the gravitons. 
As evident from the above expressions, the two-point function has no angular dependence (an outcome that follows 
from momentum conservation), whereas the higher-point functions do depend on the angles, with distinct signatures 
for the Einstein and Gauss-Bonnet corrected parts. 

This distinction between the Einstein and Gauss-Bonnet angular dependence can be attributed to the two extra 
derivatives available in the Gauss-Bonnet part of the Lagrangian. The point is that both Einstein and Gauss-Bonnet 
gravity have a two-derivative (linearized) field equation but only Einstein gravity is t ruly a two-derivative theory. 
This is what restricts and simplifies the higher-point functions of Einstein gravity [HI, |29| . On the other hand, the 
four- and higher-point functions of any Lovelock theory will generally be sensitive to the additional derivatives and, 
therefore, not subject to the same degree of constraint. Still, one can constrain any Lovelock theory to some degree 
by limiting the order of e, just like we are doing here. 



A. The gauge-theory perspective 

We are interested in the stress-tensor multi-point correlation functions for the boundary gauge theory [jig . The 
gauge-gravity duality implies that the stress-tensor correlation functions are related to the connected graviton n-point 
functions. To make the connection explicit, an appropriate process of holographic renormalization has to be applied. 
This process is straightforward in our case because, in the high-momentum regime, the connected and IPI graviton 
amplitudes agree. This will be shown explicitly for several cases, after which a general argument will be provided. 

To put this on a more formal level, we need to apply the standard rules of holographic renormalization [30-32] to the 
bulk multi-point functions. This is essentially a three-step procedure: The first step (which has already been carried 
out) is to extrapolate the bulk quantity to the boundary. [43j The second step is to multiply the extrapolated functions 
by a factor f2 9 , where is an appropriate conformal factor and the power q is determined by the conformal dimension 
of the operators whose correlation function is being calculated. The third step is to subtract off any divergences by 
using suitable boundary counter-terms (these essentially compensate for contributions from the background geometry). 
One then only retains the part that survives in the r — > oo limit. 

Let us elaborate on the second step. The conformal factor can be deduced from the asymptotic form of the AdS 
metric (see below Eq. (J7J). From this, the appropriate conformal factor can now be identified as = r/L , and one 
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then multiplies by fl q such that q = A — 3 . Here, each operator of (mass) conformal dimension A, contributes this 
amount to A = A.j, while the subtraction of 3 is meant to "strip off" the metric determinant. The operators in 

i 

our case are the gravitons for which A^ = 2 (this can be deduced from the boundary behavior of the metric) and 
the derivatives for which Ay = 1 • 

Following the described procedure, we find that all powers of r and L are stripped away from our previous expres- 
sions. For instance, 
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and similarly for the others. 

What is left is a quantity that is finite and well defined at the boundary. Hence, the first and third steps turn out 
to be of no direct consequence in the specific case that we consider. Subleading contributions to the metric could be 
important in principle but, in our case, only show up in terms that are asymptotically vanishing. 

One might still wonder about hidden implications from the third step, as the subtraction process is the only subtle 
element of the bulk-boundary dictionary. However, there are none in the high-momentum regime. To understand 
why, recall that the subtraction is tantamount to a process of matching and stripping off the (divergent) bulk and 
boundary conformal factors, and then eliminating contributions from the AdS background geometr y. T hat this is the 
underlying premise is made clear in some of the precursory works to holographic renormalization (3, I33I - I35T ] . And such 
a process would not have any bearing on our basic forms because the metric components g u and g zz are dispersed 
democratically in all our amplitudes and exhibit the same radial dependence at the boundary. Alternatively, from the 
counter-term perspective, we will see later that contributions from these have no opportunity to enter into the formal 
calculations of the graviton amplitudes. All of this simplicity hinges heavily upon the high-momentum regime being 
in effect. 



IV. A REMINDER ON COMBINATORICS OF FEYNMAN DIAGRAMS 



We include this discussion for completeness. The reader who is well versed in Feynman diagrams is advised to skip 
ahead to Eq. ([2"2"|). 

Let us start by considering a generic n-point function for a general theory. Working out the Feynman combinatorics 
amounts to (i) collecting the various terms that contribute to a Gaussian integral of the schematic form 



(17) 



and then (it) enumerating the different ways of pairing up the £'s. 

The monomial Q n represents the external legs and the exponent is an expansion of the Lagrangian in powers of (,. 
The expansion coefficients a, b, c, . . . are theory-dependent numbers that play no role in the combinatoric part of the 
analysis. The symmetrization factors are, however, important. 

Expanding out the exponent, we have 
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with the once-exponentiated monomials now representing interaction vertices. A term in this expansion makes a 
contribution to a full n-point function (C 1 )fuI whenever it contains an even number of C's, and the weight of this 
contribution is determined by counting the number of distinct ways of contracting the available £'s in pairs. In short, 
one is counting the associated number of Feynman diagrams. 

Connected and 1PI functions are more constrained than the full functions. A connected function (( n )con requires 
that two £'s from the same monomial cannot be contracted together, whereas a 1PI function (C")ip/ (or a connected 
function with no internal lines) further requires that any £ from a once-exponentiated monomial can only be contracted 
with a £ from the external monomial £™. 
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FIG. 1: Four-point function. The connected function on the left is equal to the 1PI function on the right. 



Since we are working at tree-level, even more constraints apply to the connected functions. First consider an 
arbitrary term in the expansion, which is a product of m + 1 monomials, 
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One can use the "conservation of ends", 21 + E = V [27|], to calculate the number of internal lines /, given the 
knowledge of the number of external lines E and the total number of lines ending on a vertex V . Here, E — n and 
V = YliiLiii ■ Some simple considerations then lead to the number of loops going as L = I + 1 — m [27]. Putting 
this all together, one finds 
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(20) 



where rm counts the number of occurrences of C*- The tree- level constraint means restricting to L = . 

Let us apply this general methodology to the current scenario, starting with the graviton four-point function in the 
high-momentum region. Since the three-point function is identically vanishing, the only tree-level contribution in this 
kinematic region comes from 7714 = 1 and m, = otherwise. Then n = V = 4 and so / = \ (4 — 4) = , which 
identifies this as a 1PI function. Hence, it follows that 
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must be trivially true. This is depicted in Fig. 1. 

The first non-trivial case is the six-point function. The tree-level possibilities are either mg = 1 or 7724 = 2 and, 
otherwise, m t = . The former is the 1PI six-point diagram and the latter represents the convolution of a pair of 
four-point functions, making it a 1PR contribution to the connected function. These identifications lead us to 
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as shown in Fig. 2. Here, the (hh) in the denominator represents an inverse propagator which "compensates" for 
the propagator from the internal line and cq is a numerical coefficient that measures the relative weight of the 1PR 
contribution. The next step is to determine cq. 
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FIG. 2: Six-point function. The connected function has a 1PI part and a reducible part. 

Applying the previously discussed rules, we find that cq — 10 . The 1PI case requires contracting ^ 6 £ 6 /6! with no 
pairs from the same monomial. This can be done in 6! ways, leading to 6!/6! = 1 . Meanwhile, the 1PR case entails 
the contraction of £ 6 £ 4 £ 4 /2!(4!) 2 . Drawing a single £ out of each of the C 4 's can be done in 4 2 ways. The drawn pair 
is contracted. Then the remainder C 3 C 6 C 3 /2! (3!) 2 , when subjected to the discussed rules, amounts to ( 6 ( 6 /2\(3l) 2 
and so 6!/2!(3!) 2 = 10 follows. 

If we were evaluating the Feynman diagrams for a theory of interacting scalars, this would be the end of the 
calculation. However, our actual interest is the Witten diagrams Q for a theory of spin-2 gravitons, and so things 
are more involved. 

The full calculation necessitates two additional steps. In the first, or "Step 1", we momentarily ignore all tensor 
indices and use holographic techniques to convolve a pair of four-point functions into an explicit six-point form. 
"Step 2" then accounts for the previously neglected tensor structure. We will discuss each of these in turn, beginning 
with the case of pure Einstein gravity and then the leading-order Gauss-Bonnet correction in the sequel. 

V. THE 1PR SIX-POINT FUNCTION... 
A. ...for Einstein gravity 

1. Step 1 

Let us start here with the position-space representation of the connected or, equivalently, 1PI four-point function. 
This function is depicted on the right hand side of Fig. 3. When integrated, this leads to the amplitude 
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where all tensor structure has been suppressed (as this will be accounted for in Step 2) and the holographic limit 
of r — > oo has been imposed. We also suppress labeling indices on the gravitons, with the understanding that all 
gravitons in a 1PI function and all external gravitons in a 1PR function are symmetrized. Take notice of the two 
derivatives. This number is fixed for Einstein gravity in the high-momentum regime. 

The current objective is to convolve a pair of these four-point functions and then see what it takes to manipulate 
this into an explicit six-point form. Now consider that each four-point function must "donate" precisely one internal 
graviton; otherwise, the six-point function would not be connected. Since each four-point function carries two deriva- 
tives, there are three distinct cases: (i) the two internal gravitons are both differentiated, (ii) both undifferentiated or 
(iii) "mixed" (one is differentiated and one is not). The six-point function and the different cases are shown in Fig. 4. 

Starting with the first of the three cases, we have 



(6) 



(i) 



lim 



dr 



dr / d 4 'x^-g(r,x)^-g(f,x) h(r,x) h(r,x) 



FIG. 3: Connected four-point function with two derivatives. The connected function is equal to the 1PI function. The 
differentiated gravitons are depicted as a double line and the undifferentiated gravitons, as a single line. 



V e /i(r,x) (V e /i(r,x) V//i(r,x)) V//i(f,x) /i(f,x) h{f,x) , 



(24) 



where the expectation value is used to represent the pair of internal gravitons and a superscript is included on the 
left-hand side to distinguish between the three cases. 

For the purpose of recasting this into a six-point form, we first apply integration by parts to obtain 



(6)« = 



lim 

r,f — »oc 



dr I d 4 x dr d i xy / —g(r,x)y / —g(f,x) h(r,x) 



x h(r,x) V e /i(r,x) (V e h(r,x) h(r,x)j V^|V//i(f,x) /i(f,x) /i(f,x) 
or, using the linearized field equation V a V a /i(r, x) = Q along with the product rule, 

(6)W = 



2 lim dr d x dr d x^f —g{r,x)\/ —g(r,x) h(r,x) 



(25) 



x /i(r,x) V e /i(r,x) (^V e /i(r,x) fc(f,x)> V f h(f,x) V^/i(f,x) h(f,x) . (26) 
Next, we interchange the gravitons (because they are symmetric), use the inverse of the product rule to write 
(6) (l) = -~ lim / dr / d A x [ df [ d 4 x v / -g(r 1 x) v / -g(f,x) 

3 r,r->oo J J J J 

x V e ft(r,x) ft(r,x) /i(r,x) <V/i(r, x) /i(f, x)^> /i(f , x) V//i(f, x) h(f , x) (27) 
and then integrate by parts a second time, giving 
(6) w = 



lim I dr I d 4 x dr d 4 xy^— g(r, x)a/— <?(r, x) h(r,x) 
r,r->aaj J J J 
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FIG. 4: Six-point function. Single lines denote undifferentiated gravitons and double lines, differentiated gravitons. The three 
diagrams (i), (ii), (in) differ in the way that the internal gravitons are contracted. 

x h(r,x) h(r,x) V e V e ^/i(r,x) h(f,x)\ Vfh(r,x) Vfh(f,x) h(r,x) . (28) 
Finally, we apply the Green's function form of the field equation, [45[ 



V a V a (h{r,x)h(f,x)) = = 6^(x-x)S(r-r) 
x 1 \f-9{r,x) 



, = = 5^(x-x)5(r-r), (29) 



to arrive at 



(6)W = \ lim I dr I d A x I dr I <PxyJ-g(r,-K) h{r,x) h{r,x) h(r,x) 



x S {4) (x-x)S{r-f) Vfh(f,x) Vfh(r,x) h{r,x) (30) 
or, after integrating over the (f, x) coordinates, 

(6)^ = - lim f dr f d A x\/ —g(r,x) h(r,x) h(r,x) h(r,x) h(r,x) 



3 r— >oo 



V//i(r,x) V / /i(r,x) . (31) 



This final result is the desired six-point function, containing the correct number of both derivatives and gravitons. 
Importantly, the process picked up a numerical factor of +2/3. 



11 



It is straightforward to apply the same basic procedure to the other two cases. Let us now suppose that both 
internal gravitons are undifferentiated as in the diagram marked (ii) in Fig. 4. We begin here with 



(6) (4l) = lim dr d 4 x df d 4 iV -g(r,x)y/ -g(r,Z) V e /i(r,x) 



x V e /i(r,x) /i(r,x) (h(r,x) h(r,5L)J V//i(f,x) V^/i(f,x) h(f,±) , 
then apply the inverse of the product rule 

(6><«> = \ r f^ oo f dr f d 4 x I dr [ d*xy/-g(r, x) v/- 5 (f, x) V e /i(r,x) 



(32) 



integrate by parts 



x V e h(r,x) /i(r,x) ^/i(r,x) ft(f,x)^ V//i(f,x) V^ft(f,x) ft(f,x) 



(6) ( " } = lim y dr / d 4 a; / df / d 4 i- v /-.g(r,x) v /-5(^x) V e /i(r,x) 



(33) 



x h(r,x) ft(r,x) V e ^(r,x) /i(f,x)\ V//i(f,x) V^/i(f,x) /i(f,x) , 
reapply the inverse of the product rule 

(6) (4l) = Jim /dr [ d A x f df f 'rf 4 J v qMv C ffP 



(34) 



x V e 

again integrate by parts 



ft(r,x) /i(r,x) ft(r,x) V e (/i(r,x) h(f,x)) V//i(f,x) V^/i(f,x) ft,(f,x) , 



(6) (l4) = i lim / dr / d 4 a; / df / d A x^J-g{r, x) ^- ff (f,x) /i(r,x) 



(35) 



call upon Eq. ([29 



x /i(r,x) /i(r,x) V e V e (/i(r,x) /i(f,x)) V//i(f,x) V//i(f,x) /i(f,x) , 



(6) ( " } = i lim / dr / d 4 x / df / d 4 iv/- 5 (r,x) /i(r,x) ft(r,x) h(r,x) 



(36) 



x (5< 4 ) (x - i)5(r - f ) V//i(f , x) yfh(f , x) ft(f, x) 
and integrate over the (f , x) coordinates to obtain 



(37) 



(6)(») = 1 lim / dr / d 4 xv/-ff(r,x) /i(r,x) /i(r,x) h(r,x) h(r,x) 

6 r-s-oo / / 



x V//i(r,x) V / h(r,x) 
This case comes with a numerical factor of +1/6. 



(38) 
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This leaves the mixed case (one internal with a derivative and one without) that is depicted in diagram (m) in 
Fig. 4. So we now consider 

(6) ( " l) = lim [ dr I d 4 x [ df [ d i x v / -g{r,x) v / -g(f,x) h(r,x) h(r,x) 

r.r—^oo If II 



x V e h(r,x) (S7 e h(r,x) h(f,x)j V f h(f,x) Vfh(f,x) h(f,x) , 
then apply the inverse of the product rule 



(6) ( " l) = ~ lim j dr j d 4 x f df f d 4 x v / -g(r,x) y / -g{f,x) 



(39) 



V e h(r,x) h(r,x) h(r,x) V e (/i(r,x) h(f,x)) V f h(f,x) V//i(f,x) h(f,x) , 



integrate by parts 



(6) (lri) = -i lim j dr j d A x j df j d 4 x v / -g{r,x)^-g(f,x) h(r,x) 



(40) 



employ Eq. 



h{r,x) h(r,x) V e V e (/i(r,x) h(f,x)) V f h{f,x) V//i(f,x) h(f,x) , 



(6) ( " 4) = -~ lim / / d 4 a; f (if j d A xyJ -g{ ri x) h(r,x) h{r 7 x) h(r,x) 



(41) 



x S (4) (x - x)S(r - f) V//i(f , x) Vfh(f, x) /i(f , x) 



and then end by integrating 



(6) (iri) = -i lim / dr [ d i xJ-g(r,x) h(r,x) h(r,x) h(r,x) h(r,x) 

3 r->oo / / 



(42) 



In this case, the numerical factor is —1/3. 



x V f h(r,x) V f h(r,x) 



2. Step 2 



(43) 



The next step is to determine the numerical factor that comes from the tensor structure. Let us hrst observe that 
a four-point function with two derivatives has the following form (up to additional irrelevant numerical factors and 
the metric determinant): 



(4) 1P/ = x abcd v e h ad v e h bc y™ rs h ps h 



qr j 



(44) 



where 



-\jabcd — _ I" ac bd , ab cd~\ 



(45) 



(46) 
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This structure follows from the fact that the gravitons are contracted in pairs. Undifferentiated gravitons originate from 
the expansion of the metric determinant or a contravariant metric as symmetric contractions, whereas differentiated 
gravitons originate from the expansion of the Riemann tensor as anti-symmetric contractions. The transverse/traceless 
gauge is also relevant to the stated form. 

Other than determining the placement of the tensors X and y, the derivatives play no role in this discussion, so 
we can subsequently ignore their presence (or, equivalently, work in momentum space). 

Like before, there is an important distinction between the choice of contracted gravitons, and we again start with 
the case that both of these are differentiated. Then the convolution of two four-points can be expressed as 

(6>w = yi u huh jk x ahad h ad (h bc h ps ) h qr x™ rs h tw h uv y tuvw . (47) 

We next make use of the well-known tensor structure of the graviton propagator in momentum space [36l [37| , 



(habhcd) = Qabcd = - [9ac9bd + ffadffbc — 9ab9cd] 



(48) 



We have presented the flat-space form only for simplicity. The sole effect of the AdS curvature is to alter the coefficient 
of the "trace" (or third) term in Q. As can be verified, this term never contributes in any of the cases, so we will 
subsequently drop it. It is worth noting that boundary counter-terms (as necessary for holographic renormalization; 
see Section 3.1) are non-dynamical and would therefore similarly only alter the trace. 
Hence, the above convolution becomes 



( 6 )W = y m h a h jk \x abcd h ad g bcps h qr xp^ 



Ji U -\)tuvw 
li tw 1 l uvJ / 



(49) 



Some straightforward tensor algebra then yields 



( 6 )( i > = ±y ijkl h a h jk 
= \y l]kl h a h 3k 

8 

= \y ijkl h u h jk 



h ad (s a s 5 d p + S a p S d s - 2g ad g ps )h qr X™ 
Kd (g aq 9 dr + 9 ar 9 dq + 2(d - 2) 5 Q V") h 



1 L -\jtUVW 

h h -\) tuv w 
1 l tw il uvJ / 



h ab h al 



(50) 



where h a a — has been used in the last line and d is meant as the dimensionality of the graviton sector. 

Notice that only the initial partners of the internal modes are involved in this process, with the rest of the gravitons 
simply "along for the ride". Also, the factor of 1/4 should be kept in mind. 

The other two cases follow along similar lines. If the two internal gravitons are both undifferentiated, then 



(6}<«> = x^ kl h a hJy abcd h ad g bcps h qr y pqr 



-X^ kl h u h jk 
2 



= 2^ huhj k 
= 1-X^ kl h u h jk 

And, when the internal gravitons are "mixed" , 



h ad (s a s 5 d p + 5 a p 6 d 
had(g aq g dr + g ar g dq ) 



v y prqs 



h h ytuvw 



' qr 



h n ,h ab 



(51) 



(6)(«o = y^ kl h a h 3k 



X abcd h ad Qbcps h qr y pqr 



4 J 

\y ijkl h a h jk 
4 J 

\y tlkl h a h ok 



K d (S a s 8 d p + S a p 5 d s - 2. 9 a V) V y p 
had(g aq g dr + g 



ar dq 



2g ad g q 



h ab h al 



(52) 



Again take note of the numerical factors; respectively, 1 and 1/2. 
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FIG. 5: Six-point function at order e resulting from the GB term. 
3. Putting it all together 

Let us first recall the relative weight factor of c§ = 10 that was obtained from Feynman combinatorics. One can 
see that half of this or a factor of 5 should be attributed to the mixed case, whereas a factor of 5/2 should be assigned 
to either case with a matched pair. 

We are finally in a position to determine the net relative weight of the 1PR contribution to the Einstein theory 
six-point function. The contribution of any one case is given by the product of its numerical factors for Steps 1 and 2 
times the above combinatoric factor. The three outcomes are then summed together. 

When the two internal gravitons are both differentiated, 

2 v I y 5 _ _5_ 

3 4*2 12 ' 

both undifferentiated, 

i v 1 ^ 5 5 
6 2 12 ' 

and when there is one of each, 

3 x 2 x o — 6 , 
for a total of ^ + ^ — | or zero ! 

B. ...for Gauss— Bonnet gravity 

1. Initial considerations 

Let us now consider the leading-order Gauss-Bonnet correction to the preceding 1PR calculation. By in- 
sisting on the high-momentum regime and working to order e, we must have either (hhhh) £ o (hh)~^ (hhhh) e i or 
(hhhh) e i {hh)~^ (hhhh) e o . Hence, there will be six derivatives in total and, since the process of contraction reduces 
this by two, the end result is Gauss-Bonnet's s 2 signature. Otherwise, the general procedure closely follows the 
previous Einstein theory calculations. 

2. Step 1 

At order e 1 , there are only two possible cases, as at least one internal graviton must be differentiated. These are 
depicted in diagrams (i) and (ii) in Fig. 5. Suppose that both internal gravitons arc differentiated, then the starting 



point is 



(6)«< e = lim [dr [ d 4 x [ dr [ d A x^ -g{r,x)^J -g{r,x) V d h(r,x) 

r,f->oo J J J J 



V d h(r,x) V e h(r,x) (V e h(r,x) Vfh(f,x)j V//i(f,x) /i(f,x) /i(f,x) . 
Proceeding just like before, we first apply the inverse of the product rule 

(6) w < £ = ~ lim { dr ( d A x [ dr [ d i x^/-g(r,x) v / -g{f,x) 



x V c 
integrate by parts 



V d /i(r,x) V d h(r,x) h{r,x) V e (/i(r,x) V^/i(f,x)) V f h(f,x) h(f,x) h{ 



(6) (i) ' e = -i lim / dr / d 4 a; / dr / d 4 x v / ~g(r,x) v / ~g(f,x) h{r,x) 

3 r.r^oaj J J J 



X 



V d /i(r,x) V d /i(r,x) V e V e (/i(r,x) V//i(f,x)) V//i(f,x) /i(f,x) /i(f,x) , 



reapply the inverse product rule 



(6) w ' e = i lim / dr j d 4 x j dr f d i xy / -g(r,x)^/-g(r,x) /i(r,x) 



x W d h{r,x) W d h(r,x) V e V e ^(r,x) /i(f,x)\ V/|V//i(f,x) /i(f,x) /i(f,x) 
again integrate by parts 



(6)«< e = - lim dr d i x / dr / d 4 i A /- ff (r,x) A /-ff(f,x) /i(r,x) 

3 r,r—^oo 



x V d h(r,x) W d h{r,x) V e V e (/i(r,x) /i(f,x)} Vfh(f,x) V f h(r,x) h(r,x) , 



incorporate the Green's function identity 

(6) w < £ = - lim [ dr [ d A x [ dr [ d i x v / -g(r,x) h{r,x) V d h{r,x) 

3 r,f->oo J J J J 

x V d h(r, x) (5 (4) (x - x)5(r - f) %h{f, x) Vfh(r, x) /i(f, x) , 
and finish by integrating over the (f , x) coordinates 

(6)W< e = - lim f dr [ d i x v / -g(r,x) h(r,x) V d h(r,x) V d h{r,x) h{r,x) 

Sr^ooJ J 

x Vfh(r,x) V f h(r,x) , 

which is the expected s 2 form. As usual, take note of the factor +2/3. 
Suppose now that only one internal graviton is differentiated, then 

(6)^^ = lim [ dr [ d*x [ df [ d i x v / -g{r,x)^/-g(r,x) V d h(r,x) 
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x V d h(r,x) V e /i(r,x) (V e /i(r,x) h(f,x)) /i(f,x) V//i(f,x) V//i(f,x) . (60) 



Applying the inverse product rule 



x V e 
integrating by parts 



(6) ( "^ e = i lim /dr / d 4 a; / df / d i i^J-g{r, x)y f -g(f,x) 



V d h(r,x) V d h(r,x) /i(r,x)l V e (/i(r,x) h(r,x)) h(r,x) V^/i(f,x) V//i(f,x) , (61) 



(6> ( * 4) < e = -^Um / dr /" d 4 a; / df / d 4 5v/-«?(r, x) V-.g(r, x) V d /i(r,x) 

x V d /i(r, x) /i(r, x) V e V e (/i(r, x) ft(f, x)\ h(r, x) V^/i(f , x) V//i(f , x) , (62) 



utilizing Eq. (|29[) 

(6)(»). e = -1 lim [ dr I d 4 x / df / d i i v / -g{r,x) h{r,x) V d h(r,x) 

3 r,r^oo J J J J 

x V d /i(r,x) <5 (4) (x-x)(5(r-f) V//i(f,x) V//i(f,x) h(r,x) , (63) 
and then integrating over f and x, we have 

(6)^' e = lim /dr [ d 4 x^/-g{r,x) h(r,x) V d /i(r,x) V d /i(r,x) /i(r,x) 

x V//i(r,x) V f h{r,x) . (64) 

This time the relevant factor is —1/3. 

5. S^ep 2 and beyond 

The results of Step 2 at order e are identical to those of the Einstein theory calculation. The salient points are that 
the order-e limit constrains all the background structure (X, y, Q) to be fixed at order e° (4(| and the additional 
differentiated gravitons are only "along for the ride" . And so there is a resulting factor of 1 /4 when the two internal 
gravitons are differentiated and a factor of 1/2 when only one is differentiated. 

As for Feynman combinatorics, we again recall the relative weight factor cq = 10 , but there is now an extra factor 
of 2 because of the two ways of placing the e. It is easy to see that this number should be distributed evenly between 
the two viable cases, leading to a common factor of 10-2/2 = 10. 

Finally, let us put it all together: 
When the two internal gravitons are both differentiated, 

3 X 1 X 10 = | , 

and when there is one of each, 
-3 x 2 x 10 = -3 ; 

so that, even at order e, the net contribution conspires to vanish. 

VI. GENERAL 1PR FUNCTIONS... 



We will next show that the previously observed cancelation persists for any diagram with a single internal line. 




(0 




FIG. 6: A reducible 2p ~ 2q diagram with 2p + 2q — 2n. This is one of the three 1PR diagrams contributing to the 2n-point 
function. 



.for Einstein gravity 



We begin with pure Einstein gravity and consider a 1PR function that is formed by convolving a single pair of 
amplitudes but is otherwise general, 



(h 2n )2p-2q 



(h 2 P) 1PI {h 2 ") 



IPI 



(hh) 



(65) 



where 2p, 2<j > 4 and 2n = 2p + 2q - 2 . Sec Fig. 6. 

It is a straigthforward exercise to generalize the steps of the previous 2n = 6 analysis. This is because the number 
of external gravitons only has bearing on the use of the product rule (and its inverse) in Step 1 and the distribution 
of the Feynman combinatoric factor. Nothing has changed regarding Step 2. 

Let us show exactly how this works, starting with Step 1 for the case of of both internal gravitons being differentiated, 



(2»>£U 



lim / dr J d 4 x J dr J d^sty' —g(r, x) y/—g(f, x) [h(r, x)] 



2p-2 



x V e /i(r,x) (Jv e /i(r,x) V//i(f,x)J) V//i(f,x)[/i(f,x)] 2<? 2 
The same manipulations as in S IVA II can be applied. We first integrate by parts 



( 2n >2p~2 9 = - i im I dr I <^ x \ df \ <£x\[-Qkr, x) \/~g(r, x) [h(r, x)] 



r,r— y oo 



I 2p-2 



(66) 



V e h(r, x) (V e h(r, x) h(f, x)) V f h(f, x) V/ ( [h(r, x)] 2q 2 



then apply the product rule 



(2, 



s (0 

'2p-2g 



= -(2q-2) lim j dr j d A x f df I <£xy/-g{r,x)y/-g{r,ic) 



(67) 



[/i(r,x)] 2p 2 V e /i(r,x) (V e %,x) fc(f,x)) Vfh(r,x) Vfh(f,x)[h(f,±)] 



~s-i2g-3 



(68) 
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next apply the inverse of the product rule 

(2n)g_ 2 , - - (1^1) , lh, \ / dr j d ' x J dr 'J ^^ffl 

x V e ( [h(r, x)] 2p - X ) ( V e %, x) fc(f , x)) V>(f , x) V^f , x) [h(r, x)] 29 ~ 3 , (69) 
integrate by parts for a second time 

x [/i(r,x)] 2p_1 V e V e ^(r,x) /i(f,x)\ V^/i(f,x) V>/i(r, x) [fc(f, x)] 2<? ~ 3 , (70) 
use the Green's function identity 

(2n 



2p-2q 



x <j( 4 )(x-x)<J(r-f) V//i(f,x) V//i(f,x)[/i(f,x)] 29 3 , (71) 
and end by integrating over the (f, x) coordinates 

x V//i(r,x) V f h(r,x) . (72) 

The resulting 1PR contribution appears to picks up a factor of |^Ey. But, on the other hand, had we reversed the 
roles played by p and q, the factor would rather be 2 ^~i ■ This is only an apparent ambiguity, as we have not yet 
accounted for the fact that the external gravitons should be symmetrized. This symmetrization can be implemented 
by weighing the choice to integrate first over r, x by 2 p+2q-2 anc ^ ^he cri °i ce to integrate first over r , x by 2 p+2q-2 > 
with the weights determined by the relative numbers of external gravitons (this rule will be further motivated below) . 
Hence, the correct factor becomes „ 2 f~ 1 » • 2g ~ 2 + „ • = p T g ~ 2 ■ 

' 2p+2q— 2 2p—l 2p+2q—2 2q—l P+q— 1 

Similarly, the case of two undifferentiated gravitons picks up a factor of 2p+2q-2 ' (2p-i)(2p-2) 2p+2q-2 ' 
tt: — t^t; — ?rr = 4 7 — : — ttt 2 — TT? — TT ■ The mixed case is a bit different than before, as the identity of the differentiated 

(2g-l)(2g-2) 4 (p+g-l)(p-l)(g-l) ' J 

internal graviton becomes important. [47| If the derivative is carried by the internal graviton "on the left" (or that from 
the 2p- function), then the factor goes according to 



2p-l 


1 


1 2<2-l 


2p-2 


1 P+9-2 


2p+2q-2 


2p-l 


' 2p+2q-2 


(2<?-l)(2g-2)_ 


2(9-1) P+q-1 



If the derivative is rather "on the right" , then — n/ 1 . , p ^ q 2 

o ' 2(p-l) p+q-1 

As mentioned, Step 2 gives the same outcomes as before (factors of 1/4, 1/2 and 1 for the three cases, respectively). 
As for Feynman combinatorics, it is a simple matter to generalize the six-point case; as in Section IIV1 Here, the 
Feynman relative weight is 2p ■ 2q ■ (2p + 2q- 2)!/(2p)!(2q)! = Apqw with w = (2n)!/(2p)!(2g)! . [H 

Next, consider that there is either a 1/p or l/q chance that any given internal graviton is differentiated and, 
respectively, a (p — l)/p and (q — l)/q chance that it is not. Meaning that the weight factor should be distributed 
in the following way: ~ • | • Apqw = Aw for two differentiated internals, • Srzi . 4pq W = 4(p — l)(q — l)w for 

two undifferentiated internals, ^ • • Apqw — 4(g — l)iu when there is one differentiated internal on the left and 
2^— • i • Apqw = A(p — l)iy when there is one differentiated internal on the right. 

Let us now combine the various factors: 
When the two internal gravitons are both differentiated, 
4^ 1 4w = w p±£z| 
P+q-l 4 p+g— 1 

when the two internal gravitons are both undifferentiated, 
1 P+g~ 2 vi v/^_iV»_ iY>,i — ...P+g- 2 



p+q-1 
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when only the one on the left is differentiated, 

2(9—1) P+g— 1 2 \ y I P+g— 1 

and when only the one on the right is differentiated, 

- V-Dp+^x x \ x 4 (P- = - w w£ i 

all of which sums to zero, irrespective of the values of p and q. 

Notice that any of the four types of diagrams makes a contribution that goes, up to ±w, as |^E| = 2 2^ 4 ■ This 

is relevant because it adds credence to our rule for symmetrizing the external gravitons. The point is that, from a 

momentum-space perspective, which gravitons are on which side is completely irrelevant. Hence, the net contributions 

(prior to the final summation) can only depend on the total number of external gravitons 2n, exactly what is found 

here. [H] 

Indeed, given that the symmetrization process must be unambiguous, reduce to the correct procedure when p = q 
and intermediary results can only depend on the sum 2p + 2q, it seems likely that our rule is the unique one. 



B. ...for Gauss— Bonnet gravity 



We now show that the same type of cancellation persists for the order-e contributions from Gauss-Bonnet gravity. 
Step 1 (and of course Step 2) works the same as for the Einstein gravity calculation, as the formalism automatically 
handles the situations when two undifferentiated internal gravitons is not viable (p or q = 2). What is now different 
is the distribution of the Feynman combinatoric factor, as discussed next. 

Let us assume for the time being that e goes with the 2p-point function. With e fixed, the relative weight factor 
is the same as it was for pure Einstein gravity, Apqw . What is different is the fraction of internal gravitons that 
are differentiated; this being 2/p on one side but 1/q on the other. The fraction that is undifferentiated is, 
respectively, (p — 2)/p and (q — l)/q. Meaning that the relative weight factor should now be redistributed according 
to | • - • Apqw — 8w for two differentiated internals, ^ Z ^ L ■ ■ Apqw = 4(p — 2){q — l)w for two undifferentiated 

internals, | • ■ 4pqw — 8(q — l)w for only one differentiated internal on the left and • | • 4pqw = 4(p — 2)w 
for only one differentiated internal on the right. 

The combined factors then go as follows: 
When the two internal gravitons are both differentiated, 
4^ x \ x 8w = 2u>2±2^2 

p+g— 1 4 P+g— 1 ' 

when the two internal gravitons are both undifferentiated, 

\ {p+q ^ P -l ){q -l) X 1 X 4(p - 2){q 1)W = , 

when only the one on the left is differentiated, 

- 9 , 1 n p t g ~? x i x 8(g - l)w = -2u>4^ 
2 (9—i) p+g— i 2 w i p+q~i 

and when only the one on the right is differentiated, 
~ 2(^-1) ltl-1 x | x 4(p- 2)w = ■ m 

The above sums to zero, as must also be the case when the 2q-point function carries the e. So that, even at order 
e, the same cancelations persist. 



VII. THE CASE AGAINST 1PR FUNCTIONS 



The results of the previous subsection are enough to tell us that the 1PI and connected four- and six-point functions 
are equal up to order e. 

What about higher-point functions? Of course, their 1PR contributions involve more complicated arrangements 
than so-far dealt with. For instance, the connected twelve-point function can be formed out of (a) three four-point 
functions plus a six-point function, (b) two six-point functions plus a four-point function and (c) two four-point 
functions plus an eight-point function, as well as (d) contributions with a single internal line. These combinations 
are shown in Fig. 7. Such calculations are exponentially more involved than those of the 2p-2q type. Yet, it can be 
argued on general grounds that any of these 1PR contributions must similarly vanish. 

First, consider that the generating functionals for the connected and 1PI amplitudes are the same up to a Legendre 
transformation. In either case, the generating functional contains an exponentiated action (the physical action for the 
connected generating functional and the effective action for the 1PI generating functional) and one can (in principle) 
read off the coefficients of the amplitudes from a Taylor expansion of the exponential. For a general theory, there is no 
reason for the connected and 1PI coefficients to be equal. Einstein gravity in the high-momentum regime is, however, 
different. General covariance and the restriction to two derivatives conspire to severely constrain the form of either 
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(a) 



(b) 





(c) 



(d) 



FIG. 7: Four examples of reducible contributions to the 12-point function. 



expansion [g, |29j. Essentially, all the relevant terms in both expansions depend on a single number, the numerical 
coefficient of the Ricci scalar. As the Legendre transformation has no bearing on this number, the two expansions 
have to be equal. 

We can clarify the above argument as follows: The generating functional for the Einstein theory connected ampli- 
tudes can be cast in the standard form (indices are suppressed) 



iW[J] = In 



V[h] 



Jf dx & [^g(aR+bA) + Jh\ 



(73) 



where 1Z is the Ricci scalar, A is the cosmological constant, J is an external source, a, b are numbers and (here only) 
x represents all five spacetime coordinates. The connected 2n-point function is then obtained from the relation 



(h(xi) ■ ■ ■ h(x 2n )) = (-i) 



2n+l 



s n w 



SJ(xi) ■ ■ ■ 8J(x 2n ) 



(74) 



j=o 



The basic idea in forming connected functions is that each variation of a source pulls down an external graviton, while 
the Taylor expansion of the Lagrangian provides for the rest of the structure, the vertices and propagators. 

Now consider that any of the connected functions of interest must contain exactly two derivatives and, therefore, 
must be linear in the Ricci scalar. Hence, each such amplitude contains exactly one factor of a^J—gR,. But this is 
all that it can contain because, in the high-momentum regime, the cosmological constant is of no consequence. The 
only opportunity for the constant to enter into the formalism is through an internal propagator but, as discussed in 
§ IV A 21 its presence never impacts upon the calculations. And so the connected 2n-point function can have only a 
single vertex, which is given by a^J—gTZ expanded to order 2n in number of gravitons. But, as discussed in Section 3, 
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this is precisely what defines the 1PI 2n-point function! To summarize, 

{h 2n ) E Con = a[^—gn] 



= (h 2n )?Pi ■ (75) 

Verifying this identity for all n would mean showing that every conceivable 1PR function vanishes, a futile task! 
Nonetheless, one such example is presented in Appendix A. 

A similar reasoning applies to Gauss-Bonnet's order-e corrections. In this case, the Lagrangian is aJZ + &A + 
e7£/^7p7£ /iI/7P ! where we have used a field redefinition [f| [3(1 H3| to transform the Gauss-Bonnet combination. Now, 
since these connected diagrams must contain exactly four derivatives and one power of e, they are necessarily linear 
in €\f— glZi_ lv ~ lp lZ liUlp . There is no room to include two derivatives from a Ricci scalar and, as before, no opportunity 
to incorporate the cosmological constant. In equation, 

= (h 2n )?£ ■ (76) 

However, let us emphasize that limiting to linear order in e is critical to the last conclusion. Had we, for instance, 
extended to order e 2 , then a connected function could be constructed out of a pair of Riemann-tensor-squared terms. 
This is already evident for the simplest case of an e 2 -order 1PR six-point function. Here, we start with the basic form 



r,r— >oo 



(6) e = lim dr d 4 x dr d 4 x\/— g(r, x)\/— g(r, x) Vj/i(r,x) 



x V d h(r,x) V e /i(r,x) {V e h(r,x) V//i(f,x)\ V//i(f,x) V c ft(f,x) V c /i(r,x) , (77) 
then integrate by parts 

(6) e2 = -2 lim ( dr { d A x [ df [ d i x v / -g(r,x)y/-g(f,x) V d h(r,x) 

x W d h{r, x) V e h(r, x) (w e h{r, x) h(f, x) J} V>/i(f, x) W c h(f, x) V* V c h(r, x) , (78) 
apply the inverse product rule and again integrate by parts 



(6) £ = - lim I dr I d 4 x df d A xJ — g(r, xW — g(f, x) V<j/i(r,x) 
3 r,¥—H 



x V d /i(r,x) h(r,x) V e V e ^(r,x) fc(f,x)> V//i(f,x) V c ft(f,x) V/ V c /i(f,x) , (79) 
employ Eq. (|29|) to integrate over the (f, x) coordinates 

(6) e2 = - lim / dr / d i x^/-g{r,x) V d h(r,x) S7 d h{r,x) h(r,x) 

x V//i(r,x) V c /i(r,x) V y V c /i(r,x) , (80) 
and, finally, use integration by parts followed by the inverse product rule to arrive at 

(6) £2 = -— lim [ dr [ d A x^J~g(r,x) V d h{r,x) V d h(r,x) V f h(r,x) 

lOr^rocJ J 

xV f h(r,x)V c h(r,x)V c h(r,x). (81) 

The Feynman relative weight and Step 2 bring in another factor of 10/4 = 5/2 , but the essential point is that 
this is only possible arrangement of derivatives — there can be no cancellation! Meanwhile, the corresponding 1PI 
function is provided by the next order in the Lovelock expansion, the six-derivative correction to Einstein gravity. 
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The general rule is that, at order e fc , the 1PI function is given by the 2(k + l)-derivative Lovelock extension, whereas 
1PR functions can be constructed out of any number of lower-order Lovelock terms (including the Einstein-theory 
term) provided that the total number of derivatives adds up to the same amount of 2(k + 1) (anything less would be 
in violation of the high- momentum regime) . All this is a consequence of Einstein gravity being the sole two-derivative 
theory of gravity and, thus, the only one with a built-in mechanism to prevent against the proliferation of derivatives 
in its connected amplitudes. 

Finally, let us readdress the issue of neglecting vector and scalar modes in the internal lines. Internal vector modes 
are inconsequential by the same reasoning that we applied to the internal tensors. This is because the trace term in 
the propagator is just as irrelevant for the vectors as it is for the tensors, from which the rest of the argument follows 
in parallel. 

But having scalar modes on an internal line is a different matter. For scalars, the trace term in the propagator is 
now relevant to Step 2 and thus allows for the cosmological constant to contribute to amplitudes. Nonetheless, we 
maintain that scalar modes are unphysical in the absence of any external source. 

To understand this, let us consider the simple case of a pair of three-point functions convolved into a four-point 
amplitude such that both internal gravitons are scalars. Up to an overall constant, the Einstein three-point function 
with a single scalar and two tensors goes as (with H used to indicate the scalar mode) 



(3}hhH = lim 



- ^S7 e h\7 e hH a a + S7 e hS7 f hH{ 
+ hhV e V e H a a - hhV e V f H{ 



(82) 



We next make the gauge choice V a H b a = aV b H a a , where a will eventually be set to a = 1. This gauge choice 
is consistent with the previous choice of radial gauge. Then, in similar fashion to previous calculations, we can use 
integration by parts and the inverse of the product rule to manipulate the above into a single term. For instance, 

(S)hhH = lim |(1 - a) I dr [ d^x^~9 h h V e V e iT Q . (83) 



4 1 

For a — 1 , the three-point function vanishes. So, (3)hhH is unphysical and, therefore, cannot play a role in 
the calculation of a physical quantity such as the connected four-point function. We could have, just as well, first 
computed the 1PR contribution to the connected four-point and then imposed the gauge; the same outcome and 
conclusion would persist. For higher-point functions, the procedure is expected to be more complicated but, without 
a source, the outcome should be similarly harmless. 

Our intention is to use the results of this paper in an AdS/CFT context and, therefore, a string-theory framework. 
Then we do need to consider the situation when a source for the scalar gravitons is present. More generally, we 
should consider a possible source for the scalars that originates from compactifications of string theory. In this event, 
such contributions can no longer be gauged away and, as already mentioned, the cosmological constant term could 
be relevant to the computation of the connected amplitudes. What saves us from this enormous complication is that 
the additional contributions are higher order in e. 

In the effective field theory of gravity that is induced by string theory, any such source for scalars comes with a 
coupling that scales at least as order e |38j . Hence, by introducing a scalar onto an internal line, one is suppressing the 
amplitude by another factor of e without increasing the number of derivatives. Hence, for the high-momentum region 
in particular, a physical scalar mode is subleading to the connected multi-point functions of the tensor gravitons. 



VIII. DISCONNECTED FUNCTIONS 



Although neglected so far, disconnected graviton functions can provide a useful check on consistency when it comes 
time to match our theoretical predictions with experiment. 

To see how this works, let us first consider the disconnected portion of the four-point function 

(hhhh) Dls = j§^(hh) 2 , (84) 

where the ratio of factorials accounts for the differing symmetrization factors. This is shown in Fig. 8a. 

As the above relation makes clear, one can use the value of the disconnected four-point function as an alternative 
means for evaluating the propagator. All that is required is the relevant combinatorics, so as to distinguish between 
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(a) 





FIG. 8: Full functions as sums of their connected and disconnected components. 



the disconnected and connected portions of the full four-point function. As per the discussion in Subsection IIV1 the 
"Feynman weight" of the disconnected part is determined by the number of ways of contracting h 2 ■ h A ■ h 2 such that 
connectivity is broken but without any loops. The answer is 4!/2! , (5lj | and so a relative weight (in comparison to 
the connected (equivalently, 1PI) four-point function of §j jr = jr • Hence, 



4! 
(2!)3 



(hh) 2 — (hhhh) Ful — (hhhh) Co 



(85) 



_ 1 [(hhhh) Fui - (hhhh) Cc 
{ ' ~ 3 (hh) 



(86) 



where the denominator on the right-hand side side now indicates division by a two-point function and not an inverse 
propagator. 

Also having utility is the disconnected portion of the six-point function, 



6' 6' 6' 6' 

(hhhhhh) Dls = -^^-^(hhhh)con(hh) + -^j^ (hh) c 



(87) 



where the first fraction in front of either term represents the Feynman weight (as per the previous paragraph) and 
the second fraction accounts for the symmetrization factors (as per Eq. ([84])). These are depicted in Fig. 8b. 

Now, suppose that the intent is to find an alternative means of deducing the connected four-point function. Then 
we should be looking at 

16! 16!^ 
H 41 - 2! (hhhKjcmihh) = (hhhhhh) Fu i ~ (hhhhhh) Con - ^ j^(hh) 



(hhhh) 



Con 



1 [(hhhhhh) Ful ~ (hhhhhh)cc 
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(hh) 



(hh) 2 



(89) 



IX. CONCLUSION 



To summarize, we have shown how to translate our previous results on the 1PI graviton amplitudes @ into their 
connected-function counterparts. This is an important ingredient if the gauge-gravity duality is to be used for 
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determining the corresponding stress-tensor correlations functions on the field-theory side. As explained in [l2j], such 
correlation functions would then provide the means for experimentally probing the gravitational dual of a strongly 
coupled fluid. 

Wc have shown by explicit calculation that a large class of the 1PR connected diagrams cancel off and have argued 
that this cancelation persists for all 1PR contributions in the high-momentum regime for both Einstein gravity and 
its leading-order correction. 

In many instances, a vanishing outcome where it was not expected can be the result of a symmetry principle. Is 
there such a principle here? We are looking at a specific kinematic region for which the radial derivatives are neglected 
and in essence are assumed to vanish. Away from the high-momentum regime, the radial derivatives will have some 
non-vanishing value which would correspond to spontaneously breaking the unknown associated symmetry In this 
sense, the radial degree of freedom could be viewed as a Goldstone mode. 

From the perspective of the boundary gauge theory, the radial coordinate represents an energy scale. Radial 
differentiation corresponds to a flow in this scale. That the flow has been rendered inert suggests a conformal fixed 
point of the gauge theory. Hence, the observed cancelations could indicate an unbroken conformal symmetry of the 
boundary theory for a finite N. 
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Appendix A: "Four-Four-Four" 



Here, it will be shown that, at least in one instance, a reducible function consisting of more than two components 
(i.e., more than one internal line) does indeed vanish. Out test case is the convolution of three four-point functions 
into an eight-point function, 

8 (hhhh) 1 p I (hhhh) 1PI (hhhh) 1PI 
{k >4 - 4 - 4 - Mihh) ' ( } 

We consider Einstein's theory and start by convolving the interior four-point function with one on the exterior. This 
is a different calculation than that of § IV Al because the interior four-point function has only two external gravitons. 
Hence, this is more akin to a 2p-2q (p ^= q) convolution. Also, it is now important to keep careful track of the various 
structures. 

This first convolution consists of eight distinct cases, which then leads to sixteen different diagrams depending on 
how the derivatives are initially arranged. We will work through one case in detail and then report the findings of the 
other seven. 

Let us thus consider the two diagrams in Fig. 9, 

( 6 *)(i.2) = lim f dr f d i x f df [ d A xyJ~g(r,K)^-g(?,yL) h(r,x) 

r,r-¥oo J J J J 

x h{r, x) V e h{r, x) (v e h(r, x) Vfh(f, x)^ Vfh*(f, x) h(r, x) h(f, x) , (A2) 

where the superscript on the left side denotes the particular diagram(s) and the meaning of asterisk on the graviton 
is the following: The interior four-point function contains two internal gravitons, one of which is to be contracted in 
the convolution of Eq. (|A2[) and one of which remains to be contracted in the convolution to follow. The former is, 
as usual, included within the expectation value, while the latter is marked with an asterisk for future reference. 

Starting with integration on the left (the (r, x) coordinates) and applying the usual manipulations, we arrive at two 
distinct terms, denoted as b\ and &2, corresponding respectively to diagrams (1) and (2) in Fig. 9: 



bi = /3% lim I df J d^x^J '— g(f, x) h(f,x) /i(f, x) hs(f,5i) 

V/ h* s (r, x) V/ft(f, x) h(f, x) , (A3) 
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(1) 






(2) 























FIG. 9: Two diagrams of 4-4-4 depicting the case discussed in detail in the text. The initial convolution is between the left-most 
and interior four-point functions. The difference between the two diagrams is in the right-most part of the diagram where the 
double lines and single lines are exchanged. 



f$2 Jim / df I d A x\J— g(f, x) h(f,x) h(f,Sc) hs(r, x) 



V/ V//is(f,x) h{f,5t) h(r, x) 



(A4) 



where /3\ and /?2 are relative weights which are determined below. 

Notice that the process of tensor contraction has, for either b\ or 62, induced a symmetric contraction between 
the surviving internal graviton and its new partner. (The end product of Step 2 is always a symmetric contraction; 
c /i § IV A 21 ) Because of the various integrations by parts, the distinction between symmetric and anti-symmetric 
contractions is no longer as simple as looking for which gravitons are differentiated. As far as the surviving internal 
graviton is concerned, this distinction is important in the sequel, and so we label both it and its new partner with a 
subscript of S or A accordingly. 

Next, bi and 62 are assigned the respective relative weights of 



6 = ^ 2 - - - - = — 

P1 5 3 4 2 6 120 ' 



(A5) 



3 11111 
5' ' 3 ' 4 ' 2 ' 6 ~ 240 



(A6) 



Let us explain how these weights are obtained, focusing on /3i (/?2 follows similarly). From left to right, | is due 



to our symmetrization rule, the next two numbers are a consequence of integrating by parts (2 from the product rule 
and I from the inverse of the product rule), j is the contribution from Step 2, ^ is the relative ratio for the case 

1 is the relative ratio for the case in which both of the 



in which the left-most internal graviton is differentiated and g 
right-side internal gravitons are differentiated. 

If we, rather, integrate over the (f, x) coordinates first, there is only one result, 



h = ft 



lim dr d 4 x\/— g(r), x) h(r,x) V e h(r, x) V e hs(r, x) 
r ~ > °° J J 



x h* s (r,yi) h(r,x.) h(r,x) 



(A7) 



with a weight of 



B = 2 - 2 I I I I = J- 
P3 5 3 4 2 6 180 



(A8) 
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FIG. 10: Eight diagrams of 4-4-4. Each of the diagrams represents one of the cases in the text, with the initial convolution being 
between the left-most and interior four-point functions. The right-most four-point function depicts the last term in Eq. (|A18|) . 
First row, from left to right, cases (a-c), second row, cases (d-f), third row, cases (g-h). Each of the diagrams has a partner 
diagram (not shown), where the double lines and single lines on the right are exchanged as in Fig. 9. 

Let us summarize this case using abridged notation and call it case (b) as represented in Fig. 10: 
case b 



(6* 



\(b) _ 



hhXJh (XJhXJh) Vh* hh 



1 

+ 180 

The other seven cases and their outcomes go as follows: 
case a 



= -!- / hhh s V/lg Vhh + ^— { hhh s Wh* s h h 
IzU J z4U J 

h Vh Vh s h* s hh . 



(6* 



\(«) 



VhVhh (hh) h* Vh Vh 



case c 



= i J VhVhh s h* s hh 
(6*) (c) = / / VhVhh (hh) hVhVh* 



case d 



\(d) 



^- f hhhh Vh A Vh* A + / VhVhhhh A h* A , 
60 ,/ 90 



VhVhh (hVh) Vh* hh 



= J hhh s Vh* s hVh - / hhhs VVh * s h h 



(A9) 



(A10) 



(All) 
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case e 



case f 



case g 



case h 



/ V// Yh A.s //; /,/,. (A 1.2) 



(6*) (e) = / / hhVh (vhh) h* VhVh 



^ j hhh s h* s \/h\7h - -L i /, v// v//s //; // // 



(6*) (/) = / / hhVh (vhh) hVhVh 



— [ hhhh Vh A Vh* A / VhVh h h A h* A , 

60 / A 90 ' ' 



(6*)(fl) = / / hhVh(vhVh)Vhhh 



(6*) {h) =11 VhVhh (hVh)Vhhh 



~ j \7h\7hhhh s h* s - ^ j 1 1 li VI i n hs V/» ; 



(A13) 



(A14) 



— j hhhVh Vh S h* s + — / hhhVhh s Yh* s 



+ i / // Yh Yh I: h s h s . ( Al.-.l 



(A16) 



Summing up the results of the eight cases, we find just a pair of terms, 

(6*) = — f hhhVh Y>h s h% - — [ hhhY>hh s Vhl . 
x 1 120 J 5 120 J s 

It is straightforward to show that the convolution of either of these with the remaining four-point function, 



(A17) 



(4*) = i J VhVhhs h* s + i J hh Yh A Vh* A , (A18) 

is vanishing. These calculations follow exactly as in § IVI Al with two caveats. First, for the purposes of Step 2, if 
an internal graviton is labeled with an S (A), it is treated as undifferentiated (differentiated). Second, relative-ratio 
factors have already been assigned. This is also true for the four-point function when expressed in the above form. 

Finally, since (h 8 )e-4 = from the considerations of Section 6, the preceding result allows us to conclude that 
(h 8 }con = (h 8 )ipi at least to order e°. 
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